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Abstract 

The number of configurations of the dynamical triangulation model 
of 4D euclidean quantum gravity appears to grow faster than expo- 
nentially with the volume, with the implication that the system would 
end up in the crumpled phase for any fixed k>2 (inverse bare Newton 
constant). However, a scaling region is not excluded if we allow K2 to 
go to infinity together with the volume. 
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1 Introduction 



The dynamical triangulation model of four dimensional euclidean quantum 
gravity has recently been the subject of interest 0,1, 1,1,1, 110]. 

The canonical partition function of the model at some fixed volume (num- 
ber of four-simplices) V may be defined as a sum over triangulations T of 
the hypersphere S 4 

Z{V,k 2 )= exp( K2 iV 2 ), (1) 

T(N 4 =V) 

where iVj is the number of simplices of dimension i in the triangulation T 
with fixed edge lengths. The iVj (i = . . . 4) satisfy three constraints which 
means that only two of them are independent. We have chosen N 2 and iV 4 as 
the independent variables. For comparison with other work we remark that 
if N is chosen instead of N 2 then the corresponding coupling constant k is 
related to our k 2 as k = 2k 2 . 

The weight exp^A^) is part of the Regge-Einstein action 

s = T7^rJ2 v 2 R A = K2(pN±-N 2 ), (2) 
lb7rGo a 

V 2 10arccos(l/4) 
K2 = 8Go' P = 2~, ' (3) 

where V 2 = (324/5) -1 / 4 is the volume of 2-simplices (triangles A) in units of 
V4 = 1. Numerical simulations have shown that the system at fixed volume 
can be in two phases. For k 2 > ^(N^) (weak bare coupling G ) the system is 
in an elongated phase with high (Ra) while for k 2 < /^(iV^) (strong coupling) 
it is in a crumpled phase with low (-Ra)- The elongated phase has a low ef- 
fective dimensionality and a large average distance while the crumpled phase 
has a high effective dimensionality and a small average distance between the 
simplices. The transition at k c 2 is found to be continuous with a susceptibility 
d 2 \nZ(N^K 2 )/dK 2 2 that grows with V. This opens the exciting possibility 
that continuum behavior may be found at the phase boundary |i],[2|,||,f|,|5|.[J. 

Recently, evidence was presented that the partition function Z(N^, k 2 ) 
grows faster than exponentially in N± at fixed k 2 . This result was interpreted 
as to cast doubt on the dynamical triangulation approach to four dimensional 
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gravity, because it implies that the grand canonical partition function 

Z(k 2 , k 4 ) = exp(/t 2 iV 2 - k 4 N 4 ) (4) 

T 

is ill defined (the parameter k 4 is related to the bare cosmological constant). 
We like to argue however that it is useful to study the local and ultraviolet 
properties of the model separately from the global and cosmological proper- 
ties. The nonexistence of Z(/c 2 , i^a) need not invalidate the possibility of a 
scaling regime at large N 4 . In this letter we report on our results for the be- 
havior of Z(N±, K2) and shall argue that a possible continuum limit involves 
sending also k 2 to infinity, i.e. Go ~~ > 0. 

2 Method 

Unfortunately, no set of local moves to simulate the canonical ensemble is 
known and probably no such set exists || . Therefore, we rewrite the partition 
function (P as 

Z(V, k 2 ) = exp(AS(V)) ^exp(ft 2 iV 2 - AS(N 4 ))5 N „ V , (5) 

r 

where the sum is now over all triangulations T with S A topology. In a 
simulation we can then generate configurations with Boltzmann weight 

exp( K2 iV 2 - AS(N 4 )) (6) 

and select only those configurations with N 4 = V . The precise form of AS* 
is unimportant. We have taken followed flj,]2| and used 

AS(N 4 ) = -k 4 N 4 - 7(iV 4 - V) 2 , (7) 

with 7 = 5- 10~ 4 . This is smaller than the value used in f|, 0, resulting 
in somewhat larger volume fluctuations. As the number of triangulations 
cannot grow faster than (5iy 4 )! ~ exp(5A r 4 In N 4 ) such a term ensures that 
the volume does not blow up in the simulation. The interesting question 
is now how Z(V, k 2 ) behaves as a function of the volume and the coupling 
constant k 2 . 
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The simulation of the system 

Z' = ]>>xp(/t 2 iV 2 - K4 iV 4 - 7 (iV 4 - Vf) (8) 

T 

= Y.^M^Z{N A ,K 2 )-K 4 N A - 1 {N,-V) 2 ) (9) 
allows the measurement of d In Zj <9iV 4 through the equation 

<(v, k 2 ) = d ±^A =KA + 2im) _ v) + o (A) , (io) 

which follows from a saddle point approximation for (iV 4 ) in this ensemble. 
At the same time this simulation gives us configurations with the weights (|^) 
to measure properties of (||). 

We will not expand here upon our method to generate simplicial com- 
plexes, but for definiteness let us mention that we allow all those complexes 
where no two d- dimensional (sub)simplices share exactly the same set of 
(d + 1) points. (In our criteria were somewhat different, our present 
criteria lead to numerical results in agreement with @,|6|, [7|,[3|,[5|].) 



3 Results 

Figure [l] shows the value of function of the volume at k 2 = 0, where 

the configurations at any particular volume all have the same weight. The 
points at the highest volumes are somewhat correlated as the autocorrelation 
time became of the order of the time between volume increments. The line 
is a fit to 

k\ = a + b\n(N 4 ), (11) 
where the fitting parameters are 

a = 0.82(1), 6 = 0.0323(9), (12) 

with x 2 — 24 at 46 d.o.f. We have also tried fitting the data to the converging 
functions ^ 

^ = ^a n (iV 4 )- n (13) 

n=0 
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Figure 1: The effective critical k c a for various volumes at k 2 = 0. 
for various small k and 

Kl = a + b(N 4 )- c . (14) 

The function (|13D simply does not yield any reasonable fit. On the other 
hand, the small converging power (|H|) fits with a = 1.7(9) and c = 0.06(8), 
which is consistent with zero. This also points to a logarithm. The large error 
in a is due to its large correlation with the power c. For fixed c the error in a 
would be much smaller. Both the qualitative picture and the value of b (|12]) 
agree with data reported in |7J], produced via completely independent code. 

The forms ([□]) and (|14"D are virtually indistinguishable for small powers. 
Thus, if it is indeed a logarithm, a converging small power can never be 
excluded and will even be favored if the low volume points have somewhat 
lower k\ than that of a pure logarithm. 

At the other side of the transition, at k 2 = 2.0, the same plot (not shown) 
is just a horizontal line for all volumes we have used, which went up to 90,000. 
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Figure 2: The effective critical ctS ct function of k 2 at various volumes. To 
expand the vertical scale we have subtracted a straight line — 2.3k 2 from the 
data. The values of ^(V) are indicated with diamonds. 

The fitting parameters for the logarithmic function (|TT|) are in this case 

a = 5.662(5), b = 1(5) x 1(T 4 , (15) 

with x 2 = 7.7 at 86 de grees of freedom, which is consistent with a constant 
(i.e. b is consistent with 0). So in this phase the number of configurations 
making an important contribution to the partition function does rise only 
exponentially with the volume. 

The values of k\ at various volumes and k 2 values in the region of the 
transition can be seen in figure |2|. We have subtracted 2.3^2 from the data 
for k\ to expand the vertical scale. The diamonds indicate the value of k 2 
at each volume. For the lowest three volumes these are the points obtained 
in ||, while those of 8000 and 16000 are our own. 
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Going from right to left the curves are close together and then break away 
from each other as their slope decreases. The break away point moves to the 
right as the the volume increases. Assuming that it is correct to extrapolate 



the pattern in this figure and eqs. fll^Jlq) to large volumes, we conclude that 
n c 4 first stays constant as the volume increases and then, depending on the 
value of k 2 , starts to diverge logarithmically. We also see that the transition 
between the crumpled and elongated phase moves to larger k 2 as the volume 
increases, staying near the break away points. 

There is a close relation between these curves and the average curvature, 
due to the relation 

d(N 2 ) _ d 2 In Z(V,k 2 ) _ dnl 
dV dn 2 dV 8k 2 [ ' 

This implies that the largest change in the slope of these curves will coincide 
with the susceptibility peak and that as that peak becomes narrower the 
bend in the k c a (k 2 ) curve becomes sharper. 



4 Discussion 

These results indicate that the number of crumpled configurations grows 
factorially with the volume, while the number of elongated configurations 
only grows exponentially. This would mean that at any fixed k 2 the crumpled 
configurations will always dominate for large enough volumes. 

From this it already follows that the k c 2 of the transition between the 
phases must diverge with the volume. In ref. a converging scenario is 
favored but a diverging one not excluded. In fact, it is mentioned || that 
leaving out data points with V = 500, 1000 leads to fit with k 2 diverging 
logarithmically with volume. The large k 2 region in fig. 2 seems indeed to 
suggest that V = 1000 is somewhat too low to see the asymptotic trend. 

It is an interesting question whether the bare free energy 

F(V> 2 ) = -ln £ exp(-S) (17) 

T(N 4 =V) 

is extensive at the phase boundary, dF(V, K 2 (V))/dV = const . From (^) we 
see that this would be the case if 

Ki(V, 4(V)) - pn 2 (V) - pV^l + <iV 2 >^P (18) 
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is independent of V for large V. Assuming a logarithmic dependence of k 2 on 
V the last two terms are constant so the question is whether the difference 
of the first two is V independent. We see from fig. 2 that this unlikely: in 
the range V = 2000 — 16000, k\ ph 2.4 k 2 (V) which rises faster with than 
p/t 2 ^ 2.098 k 2 . 

Expecting extensivity with the bare Einstein-Regge action is question- 
able. It consists of two very different terms oc N 2 and oc N± which will 
require different renormalization, such that the particular linear combina- 
tion k 2 (N 2 — pN^) has to be modified. In other words, the curvature term 
will mix with the volume N4 under renormalization. A natural candidate 
is k 2 N 2 — k\(N±, k 2 )N^ + const N±, where the const may be fixed by some 
normalization condition. 

We are therefore not disturbed by the fact that the average bare curvature 
is rather different from zero at the transition [|3| . One may also contemplate 
obtaining a physical curvature in terms of a physically defined metric, but 
there are other perhaps more easily accessible quantities. Elsewhere || we 
have proposed a way of measuring the renormalized Newton constant G in 
terms of the binding energy of test particles. It is however not even clear 
to us that V should be interpreted as the physical volume as measured in 
terms of a physical metric at physical scales, because V is sensitive to the 
proliferation of baby universes JK| . 

We have presented results suggesting that the bare gravitational coupling 
Go has to go to zero in a possible scaling limit. It is interesting that also 
in matrix models of 2D gravity with unrestricted topology the number of 
configurations rises factorially with the volume and a sensible continuum 
limit can only be taken by letting the bare Go go to zero in the so called 
double scaling limit Jll], [12|, |l3j . Perhaps 2D gravity in its dynamical tri- 
angulation formulation (without restriction on topology) has also a strong 
coupling phase. 

Similarly, the scenario with Go going to zero might also be applicable to 
the case of four- dimensional simplicial gravity with unrestricted topology. In 
this model it is certain that there is no exponential bound on the number of 
configurations. 
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Note 



Just as we finished this paper we received a preprint [14] presenting data that 
favors the scenario with k\ converging according to (|14]) with an exponent c 
of 1/4. 
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